Oxygen-covered tungsten crystal shape: time effects, equilibrium, surface energy and 

the edge-rounding temperature 
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The equilibrium crystal shape (ECS) of oxygen-covered tungsten micricrystal is studied as a func- 
tion of temperature. The specially designed ultrafast crystal quenching setup with the cooling rate 
of 6000 K/s allows to draw conclusions about ECS at high temperatures. The edge-rounding transi- 
tion is shown to occur between 1300 K and 1430 K. The ratio of surface free energies 7(111)/7(211) 
is determined as a function of temperature. 

PACS numbers: 68.35.Md 65.40.gp 68.60.Dv 68.37.Vj 



I. INTRODUCTION 

The study of equilibrium crystal shape (ECS) yields 
information on temperature-dependent surface free en- 
ergy, which is one of the fundamental quantities in sur- 
face physics and chemistry-. Understanding ECS is an 
important step in the ongoing challenge to construct 
realistic models of solids. The problem of ECS for 
pure one-component crystals is already complex (see eg. 
reviews^"—) . Even more challenging is the considera- 
tion of ECS with adsorption, which is important for 
the understanding of fundamental processes on catalyst 
surfaces^. There exist only few models and experimental 
reports concerning ECS with adsorptior^"— . 

The general property of ECS is that it evolves with 
temperature. At absolute zero polyhedral shapes are 
predicted2.dij with atomically sharp vortices and atomi- 
cally sharp edges. As the temperature is increased, en- 
tropic effects cause rough regions to appear on the crystal 
surface. At a certain temperature the vortices and the 
edges become rounded. We will denote by the ver- 
tex rounding temperature^ and by Tg the edge rounding 
temperature. The dependence of ECS on temperature is 
illustrated in Fig. [U (a)-(f). 

Given the difficulty of experiments on ECS with ad- 
sorption, it is reasonable to choose an exprimentally sim- 
ple and well studied adsorption system - O /W (for a re- 
view on 0/W, seei^). We treat 0/W as a model system, 
which gives insight into general properties of thermal evo- 
lution of ECS predicted by theory. On the other hand, 
0/W is one of the candidate adsorption systems for the 
fabrication of atomically sharp tungsten electrodes ap- 
plied as practical electron and ion point sourcesi^"— , and 
understanding of the temperature dependence of 0/W 
ECS allows one to control the sharpness of the electrode. 

The adsorption of oxygen on tungsten is qualitatively 
different for low and high oxygen exposures. Thermal 
desorption spectroscopy studies reveal that after low ex- 
posures (< 2 langmuir) adsorbate is all desorbed as O 
atoms, while with higher exposures oxygen is removed as 
tungsten oxides-^— . In the present work we concentrate on 
the low exposure case: 1.4 langmuir, where no alloying 
occurs. The exposure of 1.4 langmuir corresponds to a 




FIG. 1: (a)-(f) The equilibrium crystal shape in the vicinity 
of the crystal vertex depends on temperature, (g)-(i) The 
actual crystal shape observed after quenching. The observed 
shape depends on the cooling rate and may differ from the 
equilibrium shape. During quenching small steplike facets 
may be formed as shown in (g). 



coverage of less or equal to 5 • 10"'^^ molecules/cm^. 

Oxygen is known to increase the anisotropy of the 
surface free energy of tungsten and thus to induce pro- 
nounced thermal faceting of the surface (see^i and refer- 
ences therein). The O/W faceting is observed in a very 
wide temperature range: it is already observed at 800 
K2^, while annealing at 1800 K still leaves enough oxygen 
on the surface for faceting to occur. Oxygen adsorbate is 
mobile above 800 K - it forms a two-dimensional lattice 
ga!«. 

Oxygen-induced faceting of tungsten observed in ex- 
periments often takes the form of steplike (hill-and- 
valley) faceting, where steps or pyramids/pits are formed 
on the surface. These non-convex shapes are steady- 
state, non-equilibrium forms. Only under special con- 
ditions it is possible to observe a convex, equilibrium 



2 



crystal form^'^® (global faceting). It is known that the 
vertex-rounding transition occurs at Ty — 970 ± 70 K^. 

Still, there remain unsolved some interesting questions 
about the ECS of 0/W crystals. 

• As the temperature is increased, the transition from 
non-equilibrium steplike crystal surface to a convex 
equilibrium shape may be observed^^. But as the 
temperature is increased further, it appears that a 
reverse transition is taking place - a single crystal 
edge is again disassembled into a non-convex group 
of parallel edges (steplike faceting)^, apparently 
violating theoretical predictions for ECS. 

• With the present state of microscopic techniques, 
if one needs at the same time high spatial resolu- 
tion ('^l nm), high sample temperatures (~1500 K) 
and good control of surface impurities, the only way 
to observe the ECS is to observe the crystal after 
quenching (rapidly cooling) the crystal to a lower 
temperature. The question is: what is the neces- 
sary cooling rate for the atomic configuration to be 
preserved? Is the cooling rate applied in previous 
work sufficient? (400-800 K/s^-^) 

• Is it possible possible to extract from microscopic 
observations any quantitative information about 
the dependence of the surface free energy on tem- 
perature? (the problems are: microscope image 
deformations and doubts whether the cooling rate 
is sufficient). 

• Is it possible to observe the edge-rounding transi- 
tion? (Not to be confused with the vertex-rounding 
transition observed before^.) 

The main purpose of the present paper is to answer the 
above questions. In this work the crystal is mounted on a 
special support allowing for extremely high cooling rates: 
6000 K/s. As will be shown, for high crystal tempera- 
tures this does make a difference. 



II. EXPERIMENTAL 
A. The sample crystal 

The tungsten crystal had the form of a 2.6-mm long 
needle, with the cone half-angle not exceeding 10°. The 
axis of the needle was oriented along the [111] direction. 
The apex of the needle was approximately hemispherical, 
with the average radius of curvature 270 nm. In this work 
we focused on the changes in topography near the apex of 
the needle - in the vicinity of the (111) pole - where the 
crystal shape is a good approximation of the equilibrium 
crystal shape. The surface of the crystal was observed 
using Field Ion Microscopy (FIM) ^^i'^° . 

The base pressure of the vacuum chamber was 3 x 10""'^° 
Torr. The experiments were carried out at fast pace to 
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FIG. 2: Sample mounting for temperature control and ultra- 
fast cooling. 

minimize residual gas contamination. The whole pro- 
cess of crystal formation, from crystal cleaning to crystal 
quenching, took about 5 minutes. 

The sample was carefully cleaned in situ immediately 
before every experiment. Cleaning by field evaporation, 
which is applied with great success in FIM studies of 
low temperature single-atom diffusion (see egj^), is not 
appropriate for faceting studies, because at high temper- 
ature impurities from the sample shank contaminate the 
field-evaporated tip apex. For this reason we cleaned 
the sample thermally - separately for each experiment. 
This consisted of annealing at sufficiently high tempera- 
ture (1950-2200 K), where impurity-induced faceting is 
no longer observed. Surface cleanliness was verified by 
performing "blank" experiments^: the crystal was an- 
nealed at 1100 K without adsorbate, and no faceting was 
observed. 

After cleaning, the tungsten crystal was cooled to 80 K 
and exposed to oxygen (1.4 ± 0.3 Langmuir). This pro- 
vided the starting point of all experiments described in 
this work. 



B. Sample mounting for ultrafast cooling 

The main idea behind the ultrafast cooling setup con- 
structed for this work is to heat the sample while main- 
taining a very high temperature gradient in the sample 
support. When the heating stops, high temperature gra- 
dient causes efficient conduction of heat from the sam- 
ple, leading to a rapid decrease of the sample temper- 
ature. The realization of this idea is shown in Fig. [21 
Four molybdenum rods (0 — 1.5 mm) are kept at low 
temperature (80 K) and act as efficient heat sink. The 
sample is mounted on a short (6 mm) tungsten wire 
{(j) — 0.1 mm), which acts both as a resistive heater for 
the sample and as a heat conductor for sample cooling. 
Additional two tungsten potential leads (0 = 0.1 mm) 
are used to measure the resistance, and thus to monitor 
the temperature^^. 

Figure [3] shows the resulting temperature distribution 
along the 6 mm long heating/cooling wire. Note that 
along the short 3 mm distance the temperature changes 
by almost 2000 K. This high temperature gradient is de- 
sirable for ultrafast sample cooling, but it complicates the 



3 



Temperature distribution 



2000 



1500 



^1000 



500 



Ultrafast sample cooling 




FIG. 3: Calculated temperature distribution along the heat- 
ing/cooling wire for sample temperatures 1070 K and 1880 K. 
The potential leads act as additional heat sinks near x = 2 
and X = A mm. Note the extremely high temperature gra- 
dient dT/dx ~ 500 K/mm necessary for subsequent ultrafast 
cooling of the sample. 
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FIG. 4: The temperature of the sample apex (black) and of 
the sample welding point (gray) during sample cooling. 



measurement of the sample temperature. The standard 
method-^ must be modified to account for temperature 
gradients, as described in the next section. 

C. Determination of the temperature distribution 
and the cooling rate 

It is necessary to determine the spatial and temporal 
distribution of temperature T(x, t) in the sample support 
and the sample itself. Consider the segment dx of a wire. 
The total rate of energy flow dP into segment dx is equal 
to cSd^dx (c — specific heat, S = wire cross-section, d 
— density, ^ = rate of change of the temperature). The 
total rate of energy flow dP results from the balance of 
the supplied electric power dPe , thermally radiated power 
dPr and the power gain resulting from the gradient of the 
conducted power dPc — —^^dx: 

cSd^—dx = dPe — dPr TT^dx. 

at ox 

The three contributions to dP can be expressed as: 
dP^ = '^dx (Joule's law), dPr = aaT^2V^dx (Stefan- 
Boltzmann law), and ^^dx = -^{—S\^)dx (Fourier's 
law) . We arrive at the following equation for the temper- 
ature T(x, t): 

The material coefficients in this equation all depend on 
temperature: electrical resistivity p, total emissivity a, 
thermal conductivity A, sectific heat c, and density d. We 
obtain this dependence from physical tables^"—. On the 
basis of tabular data, we approximate the dependence of 
each parameter on temperature by a function of the form 
/(T) = E^-A/ O'lT^ where M, N are integers < 5. Only 



density is assumed constant over the entire temperature 
range. Note also that the wire cross-section 5 is a func- 
tion of x for the needle-shaped sample. 

The wire configuration shown in Fig. [5] is quasi one- 
dimensional, so we use the above equation, additionally 
assuming energy conservation and temperature continu- 
ity at the wire junctions. We solve the equation for 
T{x, t) by the finite-element method. An example sta- 
tionary solution is shown in Fig. [3l To determine the 
crystal temperature, we measure the resistance of the 
segment of the heating wire between the potential leads. 
On the other hand, this resistance is equal to / ^^"^J^-*-* dx 
and this value, on the basis of the full solution for T, can 
be correlated with the temperature of the apex of the 
sample. 

To determine the sample cooling rate, we calculate the 
non-stationary solution T — T{x,t). An example result 
is shown in Fig. |4l Here the heating current was turned 
off in a short time (~ 1 ms). To lower the cooling rate, 
we program the power supply to reduce the current on a 
longer time scale. 



III. RESULTS 

A. Time effects — the significance of the cooling 
rate 

In this section we demonstrate that surface topography 
observed after the crystal is quenched may depend on the 
cooling rate. The non-equilibrium structures formed dur- 
ing quenching should not be misteken for the equilibrium 
shape (see Fig. [T|) . Both the general theory of equilibrium 
crystal shape^'^, as well as the solid-on-solid model of a 
crystal with adsorbatei^ predict that as the temperature 
of the crystal is increased, the vortices and the edges be- 
come rounded. Given the sufficiently high cooling rate. 
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this rounding should be observed in the experiment. In- 
deed it has been shown that for the 0/W crystal, the 
cooling rate in the range of 400-800 K/s is sufhcient to 
observe the vertex rounding^. However, this cooling rate 
was not sufficient to observe the edge rounding. 

The special sample setup described in Sect. [IT] allows 
us to explore the influence of the cooling rate on the ob- 
served crystal shapes. We start with a simple experiment 
of annealing the crystal at high temperature: 1800 K 
(similar results are also obtained for 1700 K). The shape 
observed after quenching depends on the cooling rate, 
as demonstrated in Fig. [S] Ultrafast cooling results in 
a steplike multiple-edge configuration corresponding to 
Fig. [T] (g). Instead of a single crystal edge, we observe 3 
groups of 3 or 4 parallel edges. The distance between the 
adjacent edges L is smaller than 15 nm. The number of 
edges in a group (n) decreases to 2 for dT/dt ^ 1000 K/s. 
For cooling rates lower than 500 K/s single crystal edges 
are observed (compare Fig. [D (h)). Finally, for cooling 
rates lower than 100 K/s atomically sharp crystal vertex 
is observed (compare Fig. [l](i)). 

In Fig. [5] the strong dependence on the cooling rate 
is evident, and none of the observed shapes is the equi- 
librium shape. It turns out that cooling of the crystal 
as fast as 6000 K/s is still not sufficient for freezing the 
high-temperature equilibrium shape. It is an interesting 
question what cooling rate would be sufficient. It is clear 
that as the cooling rate increases, the number of steplike 
facets increases (Fig. [5] (b)), and the size of the facets 
decreases (Fig. [5] (c)). The equilibrium shape would cor- 
respond to L = 0. It is difficult to extrapolate from only 
four experimental points, but it is probable that the nec- 
essary cooling rate is from one to three orders of magni- 
tude higher than the maximum rate available within the 
present experimental setup. 



B. Equilibrium crystal shape 



(a) 



As shown in Sect. lIII 51 the atomic mobility of the sub- 
strate at 1800 K is so high that ECS is not preserved dur- 
ing quenching: small facets (^10 nm) are quickly formed 
(~100 ms), so it is not possible at this temperature to 
observe ECS directly. However, for lower temperatures 
surface diffusivity of the substrate is moderate and ECS 
can be observed. But if the temperature is too low, the 
opposite problem emerges: material mobility is insuffi- 
cient to overcame the kinetic barriers. As demonstrated 
previouslySi^S^ it is best to proceed as follows: first the 
crystal should be annealed at 1500 K, where hill-and- 
valley faceting no longer occurs and a globally faceted 
shape emerges, and then the crystal should be annealed 
at the desired equilibration temperature (80 s is sufficient 
for a crystal of r = 270 nm for T > 1000 K). 

What happens when the cooling rate 400-800 K/s ap- 
plied previously^ is increased by an order of magnitude? 
The results shown in Fig. [5] demonstrate that this makes 
no difference in the temperature range 1000-1300 K. The 
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FIG. 5: Different steplike-faceted structures observed after 
annealing the crystaf at a high temperature (1800 K) and 
cooling at different cooling rates. 



results for T > 1300 K are not conclusive: on the aver- 
age, ultrafast cooling leads to a slightly higher ji?!?!, as 
could be expected, but the effect is of the order of the 
experimental error. The main factors contributing to the 
fluctuations of the measured \BD\ are fluctuations of the 
residual contamination of the crystal surface, and fluctu- 
ations of the oxygen coverage. 

An important property of ECS is that it does not de- 
pend on the sample history, but only on the final anneal- 
ing temperature (reversibility). To verify this, we per- 
formed half of the experiments by approaching the equi- 
libration temperature from lower temperatures (800 K), 
an the other half - from higher temperatures (1500 K). 
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FIG. 6: The distance between the edge end points \BD\ 
(see Fig. [7]), corrected for microscope image distortions"'^, ex- 
pressed as percent of crystal radius r. Squares denote previous 
dataS obtained for dT/dt = 400-800 K/s. Diamonds and tri- 
angles denote current results for dT/dt = 6000 K/s, with the 
equilibration annealing temperature approached from below 
(diamonds) and from above (triangles). Five points deviate 
from the main trend problably due to contamination of the 
sample surface and will be discarded in energy calculations. 




FIG. 7: (a) Crystal surface in the vicinity of the [111] crystal 
pole forms a three-sided truncated pyramid, (b) Crystal cross- 
section assumed for surface free energy calcultations. 



In both cases the starting point was the globally faceted 
crystal as described above. The result is shown in Fig. [51 
Within the experimental error, the two procedures lead 
to the same crystal shape - the ECS. 



C. Calculation of the surface free energy 

As the temperature is increased, the (111) region in- 
creases at the expense of the (211) regions. According 
to the theory of the equilibrium crystal shape'^''^^, this 
means that the surface free energy of the (111) surface 
(7111) decreases relative to the surface free energy of the 
(211) surface (7211). The experimental data allow to cal- 
culate the ratio of surface free energies 7111/7211 as a 
function of crystal temperature. In order to calculate 



this ratio, we assume that the crystal has a simple, per- 
fect form depicted in Fig. [T] According to the Wulff 
theorem we have 7111/7211 — /iiii//i2ii. In the following 
paragraphs we describe how the ratio /iiii/ft.2ii can be 
experimentally determined. 

First we determine the average crystal radius r 
from the formula well known in field electron/ion 
microscopy22i2S: E = V/kr, where E is the "best im- 
age field" (22 V/nm for argon imaging), V the "best 
image voltage" , r the crystal radius, and fc w 5 (for 
microscope samples annealed at temperatures exceeding 
2000 K). The calculated crystal radius is 270 nm. 

In the next step, we analyse the FIM image of the 
crystal annealed at temperature T and determine 
the distance between the edge end points. This would be 
easy if the microscope maginfication was uniform - but it 
is nolj^, and for faceted crystals local magnification may 
vary by as much as an order magnitudc'^^. We use the 
results of very complex and time consuming ion trajec- 
tory calculations to calculate the local enhancement of 
the microscope magnification near the crystal verte»22. 
The local magnification factor turns out to be 1.14-3.26 
for crystal shapes formed at temperatures 1000-1500 K. 
The calculated distances \BD\ fall in the range 2.8 nm- 
38 nm. 

Once r and \BD\ are known, the ratio /iiii//i2ii can 
be determined from geometry - see Fig. [T) Up to the 
vertex rounding temperature, the crystal forms a sharp 
three-sided pyramid, and the distances \BD\ and \BC\ 
are zero. In this special case /iiii/ft.2ii = 1/cosa, where 
a is the angle between [111] and [211] directions: a = 
arccos(2^). For higher crystal temperatures the pyra- 
mid is truncated and \BC\ ^ ^\BD\. We assume that 
/i2ii remains constant and equal to the average crystal 
radius r, while hm decreases from the maximum value 
of r/cosa to a value of r/ cosa — |i3C|/(cotQ; -I- cot/3), 
where f3 denotes the angle between the segments AB and 
BC: f3 = n/2 — arccos(-^^). Consequently, we obtain 
the following formula for the ration of the surface free 
energies: 

7111 ^111 r/ cos a — |i3C|/(cot a -I- cot /3) 
7211 h2ii r 

which simplifies to 

Tin ^ ^ _ V3 \BD\ 
7211 2V2 12V^ r ' 

where |-BI?| is the distance between the edge end points 
corrected for the ion trajectory deflections, and r is the 
average crystal radius. 

The resulting ratio of surface free energies 7111/7211 
as a function of crystal temperature is plotted in Fig. H 
When the ratio is 1/cosa = = 1.061, the crystal 
is at the vertex rounding temperature (Tg — 960 K). 
The ratio decreases to 1.052 at 1300 K. This means 
that the anisotropy of the surface free energy decreases 
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FIG. 8: The ratio of surface free energies 7111/7211 as deter- 
mined from experiment. 
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FIG. 9: A crystal with single sharp edges is annealed at tem- 
perature T for time t = 10 s (circles) or t = 80 s (squares). 
Annealing above 1430 K destroys the crystal edge. 



with increasing temperature, in accord with the theory 

of ECS2'3. 



IV. DISCUSSION 



A. The significance of the cooling rate 



D. The edge-rounding temperature 



If we extrapolate the results shown in Fig. [5] to infinite 
cooling rate, we obtain an infinite number of infinitely 
small steps - that is, an equilibrium configuration with a 
rounded edge. This means that the crystal annealed at 
1800 K (Sect, nil A|) is above the edge-rounding transi- 
tion: 1800 K > Te. Even with the present ultrafast cool- 
ing setup, it is not possible to observe the edge-rounding 
transition directly ~ during quenching from high temper- 
atures the smooth edge reconstructs into steplike facets 
(as predicted theoretically - see in Figs. 7 and 8 of 
Ref.-iS.). However, with the present setup it is possible to 
estimate the edge rounding temperature Tg in the follow- 
ing way. First, the crystal is annealed at 1700 K (80 s), 
then the temperature is slowly lowered to 80 K (during 
~100 s) to allow for the formation of single sharp edges. 
The crystal with sharp edges is then annealed at a tem- 
perature T, to determine whether this temperature will 
destroy the sharp edge. Annealing below Te preserves 
the sharp edge. Annealing slightly above results in 
edge rounding, but the single edge may be rebuilt dur- 
ing quenching. Annealing at higher temperatures results 
in extensive edge rounding and steplike faceting of the 
rounded region during quenching. The results are shown 
in Fig. [S] From these results we can estimate the edge- 
rounding temperature from above: < 1430 K. The 
lower bound for Te is 1300 K on basis of the results of 
Sect. ImBl (up to 1300 K freezing is fast enough) . In 
conclusion, 1300 K < Te < 1430 K. 



The concept of surface free energy is most often called 
upon in two contexts: the equilibirum crystal shape and 
the thermal faceting of surfaces. It is best when the ex- 
periments are performed by observing the crystal at the 
equilibriation temperature. As mentioned in Sect. U for 
technical reasons this is not always possible, and great 
care must be taken to ensure that no significant recon- 
struction occurs when the crystal is quenched for obser- 
vation. In Sect, mil we demonstrated that the possible 
effect of the finite cooling rate can be examined by vary- 
ing the cooling rate. 

In many of the thermal faceting experiments reported 
in the literature, the observation of the surface is per- 
formed after cooling of the crystals, and the cooling rate 
is not given. At the same time the experimental results 
are often discussed in terms of zero-temperature models. 
By coincidence, this approach is often satisfactiory. The 
crystal shape formed during cooling at moderate cool- 
ing rates is close to the crystal form characteristic for 
the temperature where massive surface diffusion stops 
(for metals near 25% of the crystal melting tempera- 
ture), and this may be close to the theoretical crystal 
form at absolute zero, with sharp edges and sharp vor- 
tices. However, it is a good experimental practice to pay 
attention to the cooling rate and to watch for its possible 
influence. For example. Song et aL'^° measured the cool- 
ing rate and showed by high-temperature LEED observa- 
tion that 60 K/s available in their setup is insufficient for 
freezing the Pd/Mo crystal form at 900 K. The cooling 
rate in our setup is two orders of magnitude higher. As 
shown in Sect. IIIIl this is sufficient for freezing the 0/W 
crystal at moderate temperatures, but still not sufficient 
for freezing the ECS at 1800 K. 

The present cooling rate could be further increased by 
shortening the heating/cooling wire (see Fig. Short- 
ening of the wire not only reduces heat capacity, but 
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also reduces thermal resistance of the wire, thus lead- 
ing to a strong increase of the quenching rate. This is 
also seen from the scaling properties of the equation for 
T{x,t) derived in Sect. Ill Cl During cooling, electric cur- 
rent is switched off (drop the first term of the equation), 
and radiation plays a smaller role than heat conduction 
(drop the second term). The equation is second order 
in X and first order in t. If Ti[x,y) is a solution, so is 
T2{x,y) — Ti{ax,o?t). This means that if the wire is a 
times shorter (making the temperature gradient a times 
larger), the cooling proceeds times faster. On the 
other hand, the influence of the wire cross section S on 
the quenching rate is much weaker, because decreasing of 
the cross section area reduces heat capacity of the wire, 
but increases its thermal resistance, and the two effects 
cancel out (if radiation is neglected, S on both sides of 
the equation cancels out). 



B. Comparison with the theories of ECS 

The temperature dependence of ECS for oxygen- 
covered tungsten crystal follows the general pattern pre- 
dicted for crystals without adsorbatei"— lii. Increasing 
the temperature makes the surface free energy and the 
ECS more isotropic. At low temperature sharp edges and 
vortices are observed, and at elevated temperature the 
vortices and the edges became rounded. For 0/W crys- 
tal the vertex- rounding transition occurs at = 960 K. 
In Sections IIII Al and IIII Dl we present evidence that the 
edge-rounding transition also occurs, and estimate its 
temperature: 1300 K < < 1430 K. The increas- 
ing isotropy of the surface free energy is quantitatively 
shown in Fig. |S] - the ratio of the surface free energies 
7(111)/7(211) decreases from 1.061 at 960 K to 1.052 at 
1300 K. 

There exist only few models and experimental reports 
concerning ECS with adsorption^"—. The solid-on- 
solid model developed by Niewieczerzal and Oleksyi^ 
predicts the temperature dependence of an adsorbate- 
covered crystal shape. Although the parameters of their 
model were adjusted for palladium-covered tungsten, 
the adsorption systems Pd/W and 0/W exhibit certain 
similarities^, and the model shows qualitative agreeme- 
nent with the present work: (a) To overcame the ki- 
netic barriers and obtain ECS, the crystal must be first 
heated to a temperature where vortices and edges be- 
came rounded (Niewieczerzal and Oleksy call this ,,the 
defaceting temperature"—), (b) At low temperature the 
ECS is polyhedral, at high temperature vortices and 
edges are smoothly rounded, (c) Finite cooling rates at 
high temperatures lead to formation of non-equilibrium 
steplike facets. 



C. Implications for faceting on flat crystals 

The ECS is related to the problem of thermodynamic 
stability of flat crystal surfaces. Herring^ showed that a 
fiat crystal surface is stable if and only if the correspond- 
ing crystal orientation is present on the ECS. If a facet of 
orientation {hkl) or a rounded vertex of orientation [hkl) 
is present on the ECS, then a flat crystal of orientation 
[hkl) is thermodynamically stable. On the other hand, a 
sharp vertex implies that the corresponding flat surface 
is unstable and steplike (hill and valley) faceting will be 
observed. As the ECS depends on temperature, the sta- 
bility of some flat surfaces is also temperature-dependent. 

Song and coworkers''*''^^ observed faceted/planar phase 
transitions on flat crystals for Pd/Mo(lll) (at T = 
850 K), 0/Mo(lll) (at T = 950 K) and Pd/W(lll) 
(~1130 K, complicated by partial desorption of Pd). Ac- 
cording to Herring's theorem these correspond to vertex- 
rounding transitions for the crystal vertex in the [111] di- 
rection. Unfortunately no experiments were carried out 
for the ECS of these adsorbate-covered crystals. Pal- 
ladium adsorption and thermal faceting was studied on 
curved tungsten surfaces-^, but it turns out that a crys- 
tal of radius 200 nm cannot be thermally equilibrated in 
the temperature range where the adsorbate remains on 
the surface - kinetic constraints prevent the crystal from 
reaching a convex, globally-faceted equilibrium shape. 

It is interesting that the reversible planar-faceted tran- 
sition is observed for 0/Mo(lll) at T = 950 Ki^, 
which coincides with the vertex-rounding temperature 
r„ = 960 K for 0/W for the [111] vertex. This is unex- 
pected, because the melting point of Mo is ^-^20% lower 
then the melting point of W. Perhaps it is not the sub- 
strate properties that govern the vertex-rounding pro- 
cess, but some transition in the two-dimensional lattice 
gas formed by the adsorbed oxygen. Studies of oxygen 
on different substrates could clarify this matter. 

According to Herring's theorem, the present results 
have following implications for oxygen-induced faceting 
of fiat tungsten crystals. At the vertex-rounding temper- 
ature Ty = 960 K, a fiat 0/W(lll) crystal is expected 
to undergo a transition from {211}-faceted to defaceted 
phase, similarly to a transition observed by Song et al. 
for 0/Mo(lll)''*'. A similar conclusion can be drawn on 
basis of the edge-rounding transition. At the edge round- 
ing temperature, the edge between (121) and (211) facets 
becomes rounded, and this means that the orientation 
(332), hich is halfway between (121) and (211), becomes 
present on the ECS. For this reason at the edge- rounding 
temperature Te (1300 K < Te < 1430 K) a fiat 0/W(332) 
crystal is expected to undergo a transition from {211}- 
faceted to defaceted phase. 



V. CONCLUSIONS 

Using an ultrafast cooling setup, we studied the equi- 
librium crystal shape of oxygen-covered (1.4 Langmuir) 
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tungsten microcrystal (average radius 270 nm) in the 
vicinity of the (111) crystal pole. The results can be 
summarised as follows: 

• Up to 1300 K, cooling rates in the range ~600- 

6000 K/s arc sufficient for the observation of the 
ECS - no observable changes occur in the crystal 
shape during cooling. On the other hand, we pre- 
sented evidence that the ECS at 1700 1800 K is 
not preserved during cooling even at 6000 K/s - 
non-equilibrium steplike (hill-and-valley) facets are 
formed. 

• The ECS undergoes vertex-rounding transition at 
Ty = 960 K and edge-rounding transition at 
1300 K < Te < 1430 K. 



• The ratio of temperature-dependent surface free en- 
ergies 7(111)/7(211) decreases from 1.061 at 960 K 
to 1.052 at 1300 K. 

• A faceting/defaceting transitions are expected to 
occur for flat 0/W crystals: at 960 K for (111) 
orientation, and between 1300 K and 1430 K for 
(332) orientation. 
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